Experimental observation of the mobility edge in a waveguide with correlated disorder 
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The tight-binding model with correlated disorder introduced by Izrailev and Krokhin [PRL 82, 
4062 (1999)] has been extended to the Kronig-Penney model. The results of the calculations have 
been compared with microwave transmission spectra through a single-mode waveguide with inserted 
correlated scatterers. All predicted bands and mobility edges have been found in the experiment, 
thus demonstrating that any wanted combination of transparent and non-transparent frequency 
intervals can be realized experimentally by introducing appropriate correlations between scatterers. 
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Starting from the pioneering paper by Anderson a 
lot of progress has been achieved in the theoretical study 
of ID tight-binding models. This model includes a wide 
range of different physical situations lying in between two 
limit cases: ideal periodic lattices where all states are ex- 
tended, and completely random lattices where any state 
is exponentially localized. Specific interest has been paid 
to the so-called pseudo-random (or deterministic aperi- 
odic) potentials which demonstrate either localization or 
derealization, depending on their parameters |§-0|. A 
widely used model is described by the Harper equation 
with the site potential V n — ecos(27rcm). For a irra- 
tional, the incommensurability of the potential gives rise 
to a localization-delocalization transition (for all states) 
when the amplitude e passes through the critical value 
e cr = 2, see e.g., Ref. |J. For fixed e the energy spectrum 
of the Harper equation exhibits the famous Hofstadter 
butterfly Q when a scans the interval [0, 1]. This rather 
exotic spectrum was recently observed experimentally @] 
by making use of the equivalence of the Harper equation 
and the wave equation in a single-mode electromagnetic 
waveguide with point-like scatterers. 

For a long time a coexistence of localized and extended 
states in the spectrum of eigenenergies of ID random po- 
tentials was considered to be impossible. However, it 
was shown in Refs. |8|^] that a discrete set of delocalized 
states appears if short-range correlations are introduced 
in a random potential. This is done by repeating twice 
each value of site potential (dimer model). Recently dis- 
crete extended states have been observed in the experi- 
ment with GaAs-AlGaAs random superlattices |10|. 

A general case of ID potential in tight-binding approxi- 
mation with arbitrary correlations was considered in Ref. 
jllj . A direct relation between the pair correlation func- 
tion and the localization length has been derived. This 
relation shows that the mobility edge does exist in ID 
geometry. A few examples of potentials with correlated 
disorder were given. All these potentials necessarily con- 
tain the long-range correlations which thus give rise to a 
continuum set of delocalized states and to mobility edge. 



In this Letter we exploit the analogy between the prop- 
agation of quantum particle and electromagnetic wave, in 
order to demonstrate experimentally the existence of the 
mobility edge in ID geometry. The experimental set-up 
is shown in Fig. |l|. It is the same as has been already used 
in the microwave realization of the Hofstadter butterfly 
jjj . From the top of a waveguide of total length of 2. 15 m 
100 micrometer screws can be turned in. By varying the 
lengths of the micrometer screws different scattering ar- 
rangements can be realized. Complete information about 
the scattering matrix of a given arrangement is obtained 
via two antennas at the ends of the waveguide using a 
Wiltron 360B network analyzer. 
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FIG. 1. Experimental set-up. All dimensions are given in 



If the screws are approximated by delta scatterers, the 
propagation of a single mode through the waveguide is 
described by the wave equation for the Kronig-Penney 
model, 

oo 

i//'(z)+Ei/)(z)= EU n ^j{z n )S(z-nd). (1) 
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Here the wave function ip is associated with electric field 
of the Ti?- mode, and the energy is given by E = k 2 , 
where k is the wavenumber. We write Eq. ([!]) in the 
discrete form Jl2[ , 

tpn+l + VVi-i = [2cos(fcd) + U n kdsm(kd)} ip n , (2) 

where ip n = ip(z n )- The potential strength U n = e + £n is 
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split into two parts, its mean value e = (U n ) and fluctu- 
ations e n . Our treatment is based on the approach ]T^ ] 
which allows one to express the quantum model ([!]) in 
terms of the classical two-dimensional Hamiltonian map, 



Pn+i = (Pn + A n X n ) COS fl - x n sm(i, 

X n+ i = (p n + A n X n ) SID. fl + X n COS fl, 



(3) 



where x n = ip n is the position and p n is the conjugate 
canonical momentum. This map describes the behavior 
of a linear oscillator subjected to linear periodic delta- 
kicks. The amplitude A n of these kicks is defined as 



A„ 



ke n sin(fcd) 
sin/x 



(4) 



and the phase shift fi between two kicks is given by the 
dispersion relation for the Kronig-Pcnney model 

2 cos fi = 2 cos(kd) + kde sm(kd), < /i < ir . (5) 

The parameter fi plays the role of the Bloch number and 
the width of the Bloch band is defined by e. In this 
approach localized quantum states correspond to trajec- 
tories which are unbounded in the classical phase space 
(p, x) when n — > oo. Contrary, extended states are rep- 
resented by bounded trajectories. 

It is convenient to introduce the action-angle vari- 
ables (r, 9) according to the standard relations, x = 
r sin 9, p = r cos 9. Then the inverse localization length 
can be defined as 
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where 



r n +i 



= V 1 + A n sin(26»„) + A 2 n sin 2 9 n . 



(6) 



(7) 



This definition coincides with the standard one l^ 1 = 
(In l^n+iAM) 0] inside the allowed energy bands p[ . 
Here the brackets stand for the average over n. 

This Hamiltonian map approach turns out to be very 
effective in the study of completely disordered potentials 
jl4j as well as potentials with correlated disorder pT| , ^3| . 
In particular, in Ref. jllj the expression for the local- 
ization length for the tight-binding model with any kind 
of correlations in the potential has been obtained. Since 
the relation (||) has the form of the tight-binding model, 
one can use the results of Ref. Then the inverse 

localization length for the Kronig-Penney model reads, 



r\E) = k 2 



\) sin 2 (kd) 



sin jjL 



ip(fi) = 1 + 2 ^ Cm cos (2/J.1 



(8) 
(9) 



where £ m = (e n+m e„)/(e 2 ) is the dimensionless binary 
correlator. 

Relation (|J) is a starting point to obtain the conditions 
under which the mobility edges exist for ID random po- 
tentials. The function ip(fi) in Eq. (^|) is defined through 
its Fourier coefficients £ m . Then for a given dependence 
<p(n) the correlators can be calculated via 



it/2 

— I if ((j.) cos(2m/i) dfi. 
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(10) 



Now the problem is reduced to the construction of po- 
tentials having given correlators. Leaving mathematical 
details for a complete publication jl5| , we give here the fi- 
nal formula for set of random potentials with appropriate 
binary correlation function |16| : 
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An = — I \J cos(2mfi)dfi. 



(11) 

(12) 



Here Z n are random numbers with mean zero and vari- 
ance one. It is easy to check that the dimensionless cor- 
relator of the site potential Eq. (Jll ) coincides with the 
Fourier coefficient @. Eqs. (g)-(|12j) give an explicit 
solution of the inverse problem, since they reconstruct 
random potentials from the dependence 1(E). The exis- 
tence of the sharp mobility edges means that the function 
l _1 (E) has only finite (small) number v of derivatives at 
the corresponding energies. Then, the Fourier coefficients 
( |l0| ) decay slowly, £ m ~ mr^ v+1 ^ for m 1. This means 
that mobility edges are due to the long-range correlations 
in the random potential. Numerical data jflj for the 
step-function dependence (y = 0) of ip([i) show that the 
localization length, indeed, reveals sharp mobility edges. 
Mobility edges obtained for self-affine potentials [0 are 
also due to the long- wavelength component of correlation 
function. 

Let us now construct a random potential which results 
in the following function ip(/J.) (we chose this function to 
be symmetrical with respect to the point fi = tt/2), 



C$, < < /i < fi 2 < tt/2 
0, (i < fix] tt/2 > fi > fi 2 ■ 



(13) 



Here Cq = n/2(fi2 — fix) is the normalization constant 
obtained from the condition £o = 1. This dependence 
exhibits four sharp (v = 0) mobility edges in the first 
allowed zone. Their positions are given by two pairs of 
roots of Eq. (0) with fi = fix and fi = fi 2 . Using Eq. (|l^ ) 
one obtains (3q = 2Co(/i2 — Mi)/ 71 " an( l 



C 

(3 m — {sin(2rn^ 2 ) — sin(2m/ii)} . 
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(14) 
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FIG. 2. Transmission through the ID random sequence 
with correlation governed by Eq. (p"3|): (a) numerical results 
for N = 10 4 scatterers, (b) microwave transmission through 
an array of N = 100 scatterers (dotted line) and average over 
five different measurements (solid line), (c) microwave trans- 
mission through an array of N = 500 scatterers obtained by 
multiplicating the transfer matrices of five individual mea- 
surements. 

Experimentally it is difficult to measure the localiza- 
tion length directly. The accessible quantity is the trans- 
mission coefficient for a finite sample. Similar to the lo- 
calization length, the transmission coefficient can be 
expressed in terms of the classical map (||) ]l3| , 



(15) 
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Here r^v and r2N are the radii of the trajectories at time 
N, starting at a radius tq = 1 and angles 0q = and 
7r/2, respectively. This geometrical interpretation of the 
transmission coefficient is very useful for understanding 
its generic properties as well as for numerical simulations. 



FIG. 3. Same as Fig. |2|but for the complementary poten- 
tial, see in the text. 

A sequence of scattering strengths {e n } of the length 
N = 10 4 was generated by calculating f3 m from Eq. ( |l4| ) 
with fii/n = 0.2 and ^2/1" = 0.4, and substituting the 
result into Eq. (111]). Fig. ||(a) shows the resulting trans- 
mission for yj (e^) = 0.1 and e = —0.1. This value for 
e was obtained from Eq. (g) by adjusting the width of 
the allowed band to the experimental data. Experimen- 
tal data are shown in Figs. ||(b,c). The mobility edges 
are clearly seen near the points kd/n = 0.38,0.57,0.76 
which are the roots of Eq. (^) with fi/n — 0.4, 0.6, 0.8 at 
the r.h.s. Transmission spectrum ends at the band edge, 
kd/n — 0.91 for fi = n. Data below kd/n — 0.2 are not 
shown because of strong absorption in the waveguide at 
low frequencies. 

We have also calculated and measured the transmis- 
sion through the potential which is complementary to 
the previos one. Namely, it is transparent in the regions 
0.2 < fj, < 0.4, 0.6 < fj, < 0.8 and untransparent oth- 
erwise. For this case coefficients (3 m have opposite sign, 
Pa = ^(Mi - M2 + §), and C 2 = „,„ The 
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results are shown in Fig. 3(a) for the same parameters 
eo, e, and N as above. In this case again we got a reason- 
able agreement between analytical, numerical, and ex- 
perimental data for the positions of the mobility edges. 

For the experiment we used a segment of the length of 
N = 500 from each of mutually complementary random 
sequences in order to create the appropriate scattering ar- 
rangements in the waveguide. The e„ were mapped into 
screw lengths by identifying the minimum e„ value with 
a length of mm and the maximum value with a length 
of 3 mm. The average screw length of 1.5 mm deter- 
mines the width of the Bloch bands (?]]. This procedure 
ignores the phase shift due to the scatterers which proba- 
bly is responsible for the negative value of e (see above). 
Five measurements were performed with each random 
sequence. In one measurement a transmission through 
a realization of hundred scatterers have been measured. 
The dotted lines in Fig. ||b and ||b show the results of a 
single measurement, the solid lines represent the results 
of averaging over all five measurements. The expected 
transmission pattern is already visible. But we can do 
even better. Since in each measurement the complete 
2x2 scattering matrix S n (n = 1, .., 5) has been obtained, 
the corresponding transfer matrix T n is available as well. 
Then the total transfer matrix T is the matrix product, 
T = ^n- Thus it is possible to study the transmis- 

sion through arbitrary long sequences of scatterers with a 
set-up containing only 100 of them. However, because of 
absorption this technique is limited in our case to a total 
of about 1000 scatterers. Figs. ||c and [3)3 show such trans- 
mission spectra. In both cases the expected transparent 
and non-transparent regions are clearly reproduced. 
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FIG. 4. Microwave transmission through a random ar- 
rangement of N = 500 scatterers obtained by multiplying the 
transfer matrices of five individual measurements. 

As a check we studied in addition the transmission 
through the uncorrelated random sequence of 500 scatter- 
ers. The result is shown in Fig. |J. Here the transmission 
is approximately 100 times smaller as for the sequences 
with correlated disorder within the transparent region. 

Finally we would like to mention that Eq. (|l|) does not 



take into account waveguide effects in the device. That 
is why one cannot expect a complete similarity in the 
behavior of the transmission coefficient obtained numer- 
ically with that measured in the experiment. Neverthe- 
less, this equation allows one to calculate a correlated 
potential which being maped to the length of scatter- 
ers, reproduces any prescribed structure of transparent 
and untransparent frequency zones. This direct mapping 
turns out to be rather succesful if only the lowest mode 
of the waveguide is excited. For higher modes one might 
expect stronger influence of the waveguide effects. 
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